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ПАРАМЕТРИЧЕСКИХ ПРЕОБРАЗОВАТЕЛЕИ В ВЫСШИХ ЗОНАХ 
НЕУСТОИЧИВОСТИ КОЛЕБАНИИ 


Рассмотрены вопросы аппроксимации характеристик нелинейных резонаторов — 
нелинейных элементов параметрических преобразователей, работающих в высших 
зонах неустойчивости колебаний. Проведен расчёт коэффициентов для ап- 
проксимаций функциями гиперболического синуса, степенно-показательной, 
дробно-рациональной, степенной функциями и полиномом пятой степени. Прове- 
Ддено сравнение результатов, полученных с помощью различных аппроксимирующих 
функций, с реальной экспериментальной кривой. 
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Введение. Для развития современной техники актуальными являются ис- 
следования параметрических взаимодействий электромагнитных волновых 
полей значительной интенсивности с нелинейной средой (газовой смесью, 
полупроводниковым кристаллом и др.). Процессы, которые происходят в 
системах типа «активная среда — электромагнитное поле», лежат в основе 
функционирования практически всех электронно-волновых приборов: от 
вакуумных и твердотельных электронных устройств до лазерных и плаз- 
менных устройств. 

На протяжении последнего десятилетия значительное внимание 
научного мира уделяется изучению и анализу процессов, происходящих в 
существенно нелинейных системах. Нелинейно-параметрический эффект 
широко используется в электронике для создания малошумящих усилите- 
лей и высокоэффективных генераторов, элементов памяти, логических 
элементов: регистров, счетчиков, дешифраторов, сумматоров и функцио- 
нальных преобразователей и др. 

Одним из предметов исследования в рамках теории нелинейных ко- 
лебаний являются закономерности параметрического резонанса в высших 
зонах неустойчивости колебаний. Работы по повышению быстродействия и 
эффективности модуляции параметрических систем идут в нескольких на- 
правлениях. Одним из малоизученных, но перспективных направлений в 
разработке электронных средств является разработка устройств на основе 
резонансных схем. 

Нелинейный резонатор представляет собой колебательную парамет- 
рическую зонную систему, работающую в высших зонах неустойчивости ко- 
лебаний. 

Применение нелинейно-параметрических зонных (пазонных) систем 
(НПС) в качестве высококачественных усилителей, генераторов, преобра- 
зователей, а в цифровой технике - триггеров, логических элементов и эле- 
ментов памяти позволит не только улучшить технические характеристики 
аппаратуры, но и повысить их надежность в условиях неблагоприятного 
воздействия окружающей среды. Достоинства таких систем обуславливают 
их использование в измерительной технике, в качестве датчиков 
НЧ-диапазона, как малошумящих функциональных устройств. Необходима 
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разработка точных методик и широкого класса экспериментальных методов 
и средств исследования. Изучение свойств и закономерностей явлений в 
НПС, работающих в высших зонах неустойчивости, весьма актуально. При- 
менение подобных исследований лежит в плоскости нелинейных парамет- 
рических резонансных явлений в наноразмерных структурах. 

Нелинейно-параметрические взаимодействия распространены в 
природе очень широко и присущи многим нелинейным средам. Колеба- 
тельные процессы в НПС удобно моделировать с помощью нелинейных 
электрических цепей, в которых легко получить, описать и исследовать все 
явления и процессы, связанные с нелинейно-параметрическими взаимо- 
действиями. 

В теоретических исследованиях и практических расчётах электриче- 
ских цепей, содержащих нелинейную индуктивность, характеристику индук- 
тивности задают с помощью какой-либо аппроксимирующей функции, которая 
выражает зависимость между мгновенным значением магнитной индукции А 


и напряжённостью / магнитного поля в сердечнике. При выборе аппрокси- 
мирующей функции исходят из следующих соображений: она должна иметь 
аналитическое выражение, поскольку его легче исследовать и легче интер- 
претировать полученные результаты. Отличие аппроксимирующей функции от 
реальной экспериментальной кривой должно быть достаточно малым, а функ- 
ция должна отвечать физике исследуемого процесса, учитывая его характер- 
ные особенности. Вместе с тем вид аппроксимирующей функции должен быть, 
по возможности, простым, чтобы не усложнять анализ громоздкостью вычис- 
лений. Эти требования противоречивы, и в реальной ситуации при выборе 
аппроксимации нелинейных характеристик системы приходится идти на ком- 
промисс. Вопросам аппроксимации и их погрешностям посвящено множество 
работ, в частности работа [1]. При определении коэффициентов аппроксими- 
рующей функции Н = (В) возникает трудность в экспериментальном полу- 


чении кривой намагничивания классическими методами [2, 3]. 

Уже в первой зоне неустойчивости НПС возможны нелинейные коле- 
бания, а в зонах выше третьей наблюдаются существенно нелинейные явле- 
ния, поскольку в них используется сильная накачка [2]. При этом энергоёмкий 
элемент резонансного контура периодически находится в режиме насыщения. 
Поэтому выбранная функция аппроксимации, которая должна адекватно опи- 
сывать вольт-амперную характеристику (ВАХ) в ши- 
роком диапазоне, играет особо важную роль. 

Чтобы определить коэффициенты аппрок- 
симации, используем ВАХ для действующих зна- 

чений токов и напряжений, прямо пропорцио- 
нальную вебер-амперной характеристике магнит- 
ного материала. Схема для нахождения ВАХ маг- 
нитного материала для действующих значений 
токов и напряжений показана на рис. 1. 
Рис.1. Схема нахождения Произведём расчёт коэффициентов для 
ВАХ магнитного материа- аппроксимаций несколькими функциями: гипер- 
ен же Е болического синуса,  степенно-показательной, 
НИЯ В дробно-рациональной, степенной, степенной с 
линейным членом и полиномом пятой степени. 
Согласие расчётной и экспериментальной кривой будем уточнять по крите- 








600 


Физико-математические науки 








рию наименьших квадратов. При выводе соотношений не учитываются ин- 
дуктивность рассеивания, активные потери в цепи и потери на гистерезис 
ввиду малости этих величин и их несущественности для данной задачи. 

Постановка задачи. Исходными данными для расчётов ЯВЛЯЮТСЯ дей- 


ствующие значения тока /„ и напряжения И’. . Для аппроксимирующей 
функции конкретного вида Н = (В) требуется минимизировать сумму 


г _ 2 
Хр -1,), (1) 


где [,,[, - расчётные и экспериментальные значения действующего тока; 


р,- весовой коэффициент А-—й точки вольт-амперной характери- 
стики; и-— число точек ВАХ. 
Аппроксимация гиперболическим синусом. Пусть нелинейная зависимость 
Н=Л (в) представлена в виде гиперболического синуса Н = а5й ВВ, где 


@, В — коэффициенты аппроксимации, подлежащие определению. Используя 
второй закон Кирхгофа и закон полного тока для индуктивного элемента, под- 
ключённого к источнику синусоидального напряжения е(г)= О и зш @Е, с 


учётом сделанных выше допущений и рис.1 получаем следующие выражения: 
ав Е 
5% — =И„ зт в =е(!), 
ИИ (2) 
ПТ = Ш, 
где 5- площадь поперечного сечения сердечника; И’ - число витков 
обмотки возбуждения; /— длина средней линии магнитного поля 
в сердечнике; И„- амплитуда напряжения возбуждения; 
@- круговая частота напряжения возбуждения; #{— мгновенное 
значение тока в контуре. 
Рассмотрим совместно уравнения (2) относительно тока. Учитывая, 
что постоянная составляющая магнитной индукции в сердечнике отсутст- 


вует, из первого уравнения получаем выражение В = -ресоз Е. 
© 
Используя аппроксимацию Н = ай ВВ ‚, находим искомый ток: 
1 и 
= а @1 |. (3) 
и 5 


В полученное выражение входят коэффициенты аппроксимации @ 
и В, которые требуется определить с помощью ВАХ для действующих зна- 


чений тока и напряжения [3]. Для этого воспользуемся известным соотно- 
шением между действующим и мгновенным значениями тока: 


(4) 





а также разложением гиперболического косинуса в ряд Фурье 
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сй(х соз в!) = Г, (х)+ 5 Г; (х)соз 2, 


К=1 
где 1, (х), РЕ (х) — модифицированные функции Бесселя. 


Применив формулу Ру (у, получим выражение для 7”: 


272 
ри - 5? ВО, с0$ Е |= 
Й 5То 


ай [= 27, [5 = возим | 


Подставив это выражение в (4) и учтя, что ре и, =И\, 
[о 











после преобразований найдём зависимость между действующими значе- 
ниями тока / и напряжения И в приведённой цепи: 


[= аа @аВИ)-1, (5) 
\2 
ее Е [91 ое 
2 | 5’ 
Для точек ВАХ нелинейной индуктивности, за исключением области 
вблизи начала координат, выполняется неравенство [, (х) >>1, где 


= 29Ви . 
Используя асимптотическое разложение, справедливое для боль- 
ших значений х, для выражения /, (х) получаем соотношение 


0@)= ы 


У2лх. 


При х > Зотносительная погрешность при замене /, (х) на 


где а = 








е 
м2лх 
не превышает 6%. С учётом описанных выше соотношений из выражения (5) 
получим 


а:а 
= ехр(аВИ (6) 
в течво) 

6 Гот оя 
с погрешностью, не превышающей значения = 


бави) 200’ 





х=2аВИ. 
Определим коэффициенты аппроксимации методом выбранных 
точек по двум значениям ВАХ: // ‚ОИ; и /[,,Ц(,. Из выражения (6) имеем 


соотношения: 
ехр(аВИ,); = екр(ави,) 


0 а 
- 1/4 павВО, \/4ларЧ, 


Разделив второе уравнение (7) на первое, находим 
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1 и 
7.2. РВ, 0). (8) 


Из выражений (7) и (8) получаем выражения для коэффициентов 


аир: 





: ,(: [0 | а=--[лаВО ев). 
(22 


а, -0,) 1 О, 


Более точные значения коэффициентов аппроксимации можно по- 
лучить из условия минимума суммы (1), которая является функцией иско- 
мых ди р и имеет вид 


Е(а, В)= рр. ав -Т)-Ь,. 
дЕ(а, В) 
да 


Из условия = 0 находим для (© как функции В: 


Хр, уро, 


а=^ : (9) 
ор, (2ави,)-Ц 





Из условия = 0 получаем уравнение 


дЕ(а,, В) 
ОВ 
ри 24ВИ, ) _ 
баро,)-1 

Подставив в уравнение (10) выражение для &, взятое из выражения 
(9), получим нелинейное уравнение для определения р: 


ов) ру (2аВи,) |2 Хр (2ави,)- 


г ры (24ВИ») ху [1 (2ави,)-|=0. 9 

ны к (2авВ0,-1) = 

Таким образом, из выражения (11) можно найти Д, а затем из уравне- 

ния (9) можно найти @ . Если используемая часть нелинейной характеристи- 
ки индуктивности далека от области насыщения, то выражения для опре- 
деления коэффициентов х и В можно существенно упростить. Обозна- 





а. У ри @чви,) ао (10) 





чим поправки к коэффициентам ©, и Ву, найденным по методу выбран- 


ных точек, через Да и ЛВ соответственно и проведём линеаризацию вы- 


ражения (5) относительно этих поправок. Тогда для каждой точки ВАХ 
можно написать условное уравнение 








ада +ь дВ=И,, (12) 
[2 
где а, =а\/1 (248 И,)-1; В, аи = -аа,. 
0 0-7 
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Искомые поправки ло и ДВ можно найти из системы нормальных 
уравнений: 


да») рьаь + дву. Рьа, В, = У рнаи,, 
к к к 


да» рьа4Ь, + АВУ` р = э. - 
К=1 К=1 К=1 


Процесс приближения можно повторить, взяв в качестве 0, и Во 


улучшенные значения коэффициентов. Этот процесс сходится, если невяз- 
ки условных уравнений (12) не будут слишком велики. 


Начало 


Ввод исходных 
данных и, №, (к, 
м1, Р,4 


Вывод исходных Вычисление ДА, Дх, 
НЫЕ Ду, х, у, а, В, бо, 6х, 


У 









































Вычисление 




































































оо и Во 
К=1 
Вычисление 
ак, Бь Ал, Ао, Вт Вычисление 
Вэ, Су, С>, 5% Лок, Пк, [к 
Да Вывод 
результатов 
КЕК+1 





























Рис. 2. Алгоритм определения коэффициентов аппроксимации 
функцией гиперболического синуса 


На рис.2 представлен алгоритм определения коэффициентов ап- 
проксимации гиперболическим синусом по И аппроксимирующим точкам 


за ЛМ циклов, где х и у - погрешности вычисления о, и В. 
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Аппроксимация степенно-показательной функцией. Аппроксимация 
степенно-показательной функцией Н =аВехр(ВВ”) проводится анало- 


гично предыдущему случаю. Находим зависимость между действующими 
значениями тока и напряжения. Мгновенное значение тока в контуре имеет 
вид 


1= И окр (ВВ } (13) 





Из выражения (13) с учётом (3) получаем соотношение для р 
Р =(а- ааи) (1+ соз2%#)ехр(Ва?И? +Ва?и? с0з265). 


При вычислении действующего значения тока по формуле (4) вос- 
пользуемся разложением показательной функции в ряд Фурье: 


ехр(Ва?И? 605201) = (Ва?0? ) + 21, (В4?0? соз2 ог 
К= 


Сделав несложные преобразования и вычислив интеграл, получим 
зависимость между действующими значениями тока и напряжения: 


а + с0$26/)ехр(В - 470? с05261] ы (8 20? Ей (в +20? } 
0 


[=@а`адИ п (в-аео?)+п(в-аео?)еч[ 3-40? (14) 


Используя асимптотическое разложение, учтём, что для больших 





значений х верно соотношение |, (х)=1 (х)= 


АУ 


° Тогда (14) пере- 
пх 


пишется в виде 


1=а:а\/аИ еж) 
пВ 


Из этого выражения по методу выбранных точек находим: 


2 2 
[1 =а-а\/4(, В ехр(Ва?0! ); раза 9ыь _— ехр(В4*И; ). 
Отсюда получаем приближённые значения для < и В: 
1 
Е: ОМИ Ь Ч ‚ а= 1 {| ехр(-4?в0?). 
Ф (0-0?) Ци уу, и \2 
Сумма (1) теперь будет иметь вид 


Рбе.в)= тр азов атом О 1 


Используя формулы для производных от модифицированных функ- 


ций Бесселя, из 0Е (В) =0и 0 (а,В) — 0 находим более точные значе- 


ба, ОВ 
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ния коэффициентов с и В — нелинейное уравнение для определения более 
точного значения В и формулу для определения с: 


1 
дети 
о Нч) (1 доз МАРИ] 
Ч Ч» 


Ур 
5/6 ва?0 + Цао?) 


ар ехр(Ва?0? )* фил | - и ОХ 














ов) уро [обито (мое увчто? | 
> 4 2712 1 2:5 о 
+: бы [1-м (Ва И; Оо 0: |- 
-{лоьбыноя = емо офи 


о РИ ечо[ 5 в?0? 
р (Ва?) + 1 (в9*0?) 
Коэффициенты условного уравнения (12) имеют вид: 
а; =а40, оч [ВыРО: Мыло: +1 (Ва? ), 1, =1, — ола. 


1 зв (1. ва\ (89970) +1 (ВО) | 
09 Итехр >Во9 и 2772 2112 7 И 
1, (Вьч 0; +1 (Вьч 0+) Ро С 


Дробно-рациональная аппроксимация. Дробно-рациональная аппрок- 








2772 1 2712 
бы сз) * [1 зажат (840; ) = 0. 








В ы .2 
и >. После преобразований для значения 1 
В 





симация имеет вид Н = 


2 2 
4(аа4( )`с0$`@! ‚ Вычислив интеграл 


получаем выражение 22 _ 
2 
(В — 074? сов? ©!) 





ть 2 
1 Р= (чад) ‚ Получим формулу ДлЯ действующих значений тока 


® 8-90) 
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ааа 
{В(В- 970?) 





и напряжения: 7 — . По методу выбранных точек получаем 


44, _ а, 


{Ве-ги")°° Вв-ео’). 


Из этих выражений получаем приближённые формулы для опреде- 


выражения для /} и /,: и = 























ления си В: 
4 
4 НВ (:} [= 0? 
Ва] и1 о |1 {= = ее 
БО, БО, Т. и ы 
ЗП -2х | —] 
1 0, 
1 1 4 2тт2 
а=—— —-а4. 
= $1 8(в-ч°0?} 
Сумма (1) В данном случае имеет вид 
п 1 23) ь 
ибо) ур аи (8-40: -, 
&= 
Из условий 92 (%В) _ | и ОЕ (о,В) —0 получаем формулу для 
да, ОВ 


определения с, и нелинейное уравнение для определения В: 
: : п РО ь (48-4°0;) 
РИ, В “(В- 402] 4 155 (в _ 12772 
2 вЫьИ ( — и ЧВ° (В РИ: ) 
1 3 2 2772 \ 
у — -> "” р.И; (4В-а (Ц 
аа Ур, О%В *(В-4°0:) зай — | ты :) 
ыы +1 В(В-а?0;) 
РИ, _ РА: (48-90) 








0 = 














5 488-90?) Я ув(в-е0?) 


Выражения для коэффициентов условных уравнений имеют вид: 


а; =а4 0, "Во (В =, у 


го 
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В, = очачй, (В, - и: |* +3] В, (в - па 


= -оа,. 


Степенная аппроксимация. Степенная аппроксимация Н =аВ?**, 

где и=1,2,3...; ие №. Выражение для квадрата мгновенного тока имеет 
т 

вид Г =2(аа)” (40 созаг)"". Для вычисления интеграла [74 вна- 


0 
2п 


чале нужно вычислить интеграл у, = [соз” хах. Интегрируя его по 
0 


21-1 г Е 
ой 2"? ОЕ бо: = [ соз хах. Ис 
ГИ 


пользуя эту рекуррентную формулу, находим выражение для У: 
1.3.5...(2и-3)(2п-—1) р (2п)! ь Г(2" +1) 
п 


=2л =иД — . 
7 2.4.6--(2п-2).2п 2?" (п!)' 2" [Г(и+1)] 


где Г(2п +1), Г(п +1) — гамма-функции, таблицы которых имеются в 


математических справочниках, например, в работе [4]. После преобразова- 
ний получаем следующую зависимость между действующими значениями 
тока и напряжения: 


п %ч [2(4и+2) 1 2+1 
та) [= а =. 

гле -2, | +1! _2| Г(а+2) |№. 
а Р@н+2! ча рГ@н+з) 


Сумма (1) в случае степенной аппроксимации имеет вид 


у и 2и-1 
Е(а)= > р, «(8 -Г, 


частям, получаем: Ут = 











= 8 
дЕ(а) 
Из условия минимума суммы =. =0 находим выражение для с: 
2и-+1 < 2и-+1 
ы. а Г, 
А. (15) 
4и-+2 

ре 


К=1 
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Степенная аппроксимация с линейным членом. Степенная аппрокси- 


. ар 
мация с линейным членом имеет вид Н = “.В“*" + ВВ . Такая аппрокси- 


мация позволяет улучшить согласие между расчётной и экспериментальной 
кривыми в линейной части ВАХ по сравнению со степенной аппроксимаци- 


ей, в которой В =0. Зависимость между действующими значениями тока / 
и напряжения (здесь можно представить в виде 


Т=а аи" а? +п(аи >" ав + (а В?, 
Г(4и +3) | Г(2. +3) 


И 
2 [Г(2.+2)] 2% [Г(и+2)] 
Это выражение можно записать иначе: 


1=а[Д аи)" а+чив) 2, где о 60. (16) 
[а (0) +В] 


Непосредственный расчёт показывает, что при любых значениях 
О и ци величина 2 намного меньше единицы. Уравнение вида (16) реша- 


ется методом последовательных приближений, причем величину 2 в пер- 
ВОМ приближении полагаем равной нулю. Если степень многочлена 


Н=аВ?"* +ВВ достаточно высока, то коэффициент В предпочтитель- 


нее определять не совместно с коэффициентом (©, а независимо, по ли- 
нейной части данной ВАХ. Тогда коэффициент © найдётся в несколько 


где 





приближений из выражения (15), где вместо / ; Следует использовать вы- 
ражение 


1 
р 
а О. 
1-2 
Аппроксимация полиномом пятой степени. Нелинейную зависимость 
Н =(В) можно аппроксимировать полиномом пятой степени: 


Н=АВ+КВ +В, где К,К.,,К.- коэффициенты аппроксимации [3]. 
После аналогичных преобразований получим выражения для действующих 
значений тока / и напряжения (: 
ты ЗА. Е 10% 
5’ © 47 (5'®) ^ 16. 
Полученное выражение представляет собой аналитическую связь ме- 
жду действующим значением напряжения (и тока / для нелинейной ин- 
дуктивности при аппроксимации нелинейной характеристики намагничивания 
Н=Х (В) полиномом пятой степени. Получение коэффициентов аппрокси- 


Г=о0 + 0 * + 950 °; ГДЕ а, = 





мации полиномом пятой степени подробно описано в нашей монографии [3]. 
Алгоритмы определения коэффициентов аппроксимации всех рас- 
смотренных функций аналогичны алгоритму, представленному на рис.2. 
Сравнение полученных аппроксимаций с экспериментом. Описан- 
ные выше аппроксимации дают лучшие (в смысле критерия точности) зна- 
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чения коэффициентов, которые можно получить путем варьирования вы- 
бранных точек, задающих область, и целиком содержащую внутри себя 
экспериментальную кривую или интересующую ее часть. На рис.3 показа- 
ны экспериментальные и расчётные кривые намагничивания, пересчитан- 
ные к току и напряжению, для рассмотренных выше аппроксимаций, рас- 
считанные с помощью написанной на языке программирования 
ОЕЁРНТ 2006 компьютерной программы АРРКОХ1МАТА, реализующей алго- 
ритмы расчёта коэффициентов всех изложенных выше аппроксимаций. 














0,01 0.02 (0,03 0,04 0,05 0.06 0,07 0,08 ГА 
) 





а) 














0,01 0,02 0,03 0,04 0,05 0,06 0,07 0,08 ГА 





6) 


Рис. 3. Кривые намагничивания, пересчитанные к току и напряжению, полученные 
с помощью программы АРРВОХМАТА: а) 1 -— по экспериментальным данным; 
2 — аппроксимация гиперболическим синусом, Н = 0,244.10“ 5й (37,58); 3 - ае- 


пенно-показательная аппроксимация, Н=1,077В ехр(30,528°); 4 — степенная 
аппроксимация, Н=6,63.10° В?'; 6) 1 - по экспериментальным данным, 


5 — дробно-рациональная аппроксимация, р — 7,95 В ; 6 — степенная аппрок- 
0,22 — В? 
симация с линейным членом, Н = 5,25. 10° В?! + 66,8В 
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Выводы. Оказывается, что кривые, соответствующие аппроксимации ги- 
перболическим синусом и степенно-показательной функцией, практически 
совпадают друг с другом, если в обоих случаях коэффициенты © и В оп- 


ределяются по одним и тем же выбранным точкам вольт-амперной харак- 
теристики. В случае глубокого насыщения материала сердечника эти кри- 
вые хорошо согласуются с экспериментальной ВАХ в нелинейной области и 
хуже - в линейной. Этот же вывод относится и к степенной аппроксимации, 
где 2и +1=21. Заметно лучшее совпадение с линейной частью экспери- 


ментальной ВАХ наблюдается у кривой, соответствующей аппроксимации 
степенным двучленом с линейной частью, степень которого равна 21, а 
также дробно-рациональной аппроксимации, что подтверждает вывод, 
сделанный относительно неё в работе [3]. 

Таким образом, именно аппроксимация гиперболическим синусом в 
случае глубокого насыщения материала сердечника лучше согласуется с 
экспериментальной ВАХ в нелинейной области. Кроме того, эта аппрокси- 
мация приводит к получению более компактных выражений, что уменьша- 
ет громоздкость расчётов и делает именно эту аппроксимацию наиболее 
пригодной для описания сильно нелинейной характеристики параметриче- 
ской системы. 

Формулы для коэффициентов вышеуказанных аппроксимаций най- 
дут применение при расчётах характеристик нелинейных элементов пара- 
метрических преобразователей, работающих в высших зонах неустойчиво- 
сти колебаний, содержащих существенно нелинейные элементы в составе 
нелинейных резонаторов различной физической природы. 
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